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Abstract 

A method is presented for the calculation of all exact ground states of diluted 
antiferromagnets and random field systems in an arbitrary range of magnetic fields 
B e [B s tart, -Bend] resp. A £ [A s t art , A end ] . It works by calculating all jump-fields B, A 
where the system changes it's ground state. For each field value all degenerated ground 
states are represented by a set of (anti-) ferromagnetic clusters and a relation between 
the clusters. So a complete description of the ground state structure of these systems is 
possible. 

Systems are investigated up to size 48 3 on the whole field-range and up to 160 3 for 
some particular fields. The behavior of order parameters is investigated, the number 
of jumps is analyzed and the degree of degeneracy as functions of size and fields is 
calculated. 

Keywords (PACS-codes): random systems (75.10.N), ground state calculation 
(75.40.H), order parameters (75.40. C), graph theory (02.10). 

1 Introduction 

The behavior of diluted Ising antiferromagnets in a homogeneous field (DAFF) and of random 
field Ising ferromagnets (RFIM) is still not well understood (for a review see §). Although 
it has been argued that DAFF and RFIM belong to the same universality classes j^, [| recent 
results || ||, ^, [?J suggest that there might be essential differences. 

In H integer optimization algorithms where applied to the RFIM in order to calculate exact 
ground states at fixed fields in polynomial time. In |^] different degenerated ground states 
of DAFFs and RFIMs were calculated by using a simple extension of these optimization 
methods. 

The invention of new fast algorithms |l0| inspired us to develop an algorithm which calculates 
the ground states for all fields and is presented in this paper. Additionally the complete 
degenerated structure of a ground state is computable at once, by the use of some more 
complex algorithms from graph theory (jll| . Applying these methods a far better insight 
into the nature of ground states of disordered systems is possible. 

All systems we investigated were cubic LxLxL lattices of Ising spins &i — ±1 with periodic 
boundary conditions. For a DAFF each lattice site is occupied (e^ = 1) with probability p 
(otherwise ti = 0). The DAFF is described by the following Hamiltonian 
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H = J ^2 ( ,t ,n,o , - B^^eiOi (1) 

<i,j> i 

with the uniform external field B > and the interaction constant J > 0. The sum < i,j > 
runs over all pairs of nearest neighbors. For low temperatures and small fields B < B c the 
DAFF has an antifcrromagnetic phase. In this paper we show for the critical field B c sa 0.9. 
For finite fields and low temperatures exists a frozen domain state which is characterized by 
fractal domains || 0, For higher fields or higher temperatures the DAFF is paramagnetic. 
The Hamiltonian for the RFIM is 

h = -j wj -J2 Bi(Xi ( 2 ) 

<i,j> i 

All sites ore occupied, but the magnetic fields Bi are site dependent and distributed accord- 
ing a bimodal (±A-RFIM) or a Gaussian (Gaussian-RFIM) probability distribution. Their 
probability density functions are 



p(B,) = i(J(Bi-A)+J(B, + A)) (4) 

Like the DAFF the 3d RFIM has a long range ordered low temperature and low field phase, 
a disordered phase for finite fields (critical field: A c = 2.35 ||) and is paramagnetic for high 
temperatures or high fields. 

The outline of the paper is as follows. The first section explains how all degenerated ground 
states of a system at a fixed field value can be represented by a relation of clusters. The 
second section gives an algorithm which calculates the ground states for all fields. In the third 
section we show the results of our calculations. We present the order parameter for sample 
systems of different sizes, estimate the values for the critical fields, analyze how the number 
of jumps increases with increasing system size and investigate the degree of degeneracy as 
function of field and size. In the last section we summarize our results. 



2 Calculation of a ground state 

We used well known algorithms from graph theory [l^] to calculate the ground state 

of a system at given field B, A and interaction constant J. All quantities B, A, J have to be 
integer- valued. The calculation works by transforming the system into a network |jq] , and 
calculating the maximum flow [l7| . [] All degenerated ground states of the system are 
given Q by a set V = {S — Vq, Vi, . . . , V n , V n+ i = T} of clusters, and a binary relation R 
defined on V. Each cluster Vi is a set of antiferromegnetically (DAFF) resp. ferromagnetically 
(RFIM) ordered spins. These spins are not necessarily spatial connected. 
From the pair (V, R) a ground state is calculated in the following way. Each cluster \% 
(k = 0, . . . , n + 1) is assigned a orientation (or sign) s(Vfc): 

The signs of the clusters S,T are determined: s(S) — 1, s(T) = — 1. For all other clusters 
T4 (fc = 1, ...,n), called inner clusters, the sign can be chosen free (s(Vk) — ±1) under 
consideration of the condition (I — 1, ... ,n): 

if {s(V k ) = +1 and (V k ,Vi) G R) then s(Vi) = +1 (5) 

I.e., if a cluster has the sign +1, then all successors in the relation R must have the sign +1 
too. 

1 Implementation details: We used Tarjan's wave algorithm together with the heuristic speed-ups of Traff. 
In the construction of the level graph we allowed not only edges (v,w) with level(io) = level (v)+l, but also 
all edges (v, t) where t is the sink. For this measure, we observed an additional speed-up of roughly factor 2 
for the systems we calculated. 
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From the signs of the clusters the spin states of the spins i £ Vk are calculated by 

<Ti= t(i)s(V k ) (DAFF) (6) 
a l= s(V k ) (RFIM) (7) 

The function t(i) divides the lattice of the DAFF into two sublattices: t(i) = (— \) x +y +z for a 
spin i with coordinates (x, y, z). So the clusters S, T consist of the spins with fixed states, for 
each direction one cluster, all other clusters contain spins which contribute to the degeneracy 
of the ground state. The dependencies between spins of adjacent clusters are given by the 
relation R 

As example we consider an one dimensional RFIM chain with open ends consisting of four 
spins 0i, . . . , 04. The random fields are B\ — 2, B 2 = B 3 = 0, B4 = — 2 (J = 1). So the 
states of the end-spins are fixed in the ground state <j\ = 1,0-4 = — 1- The states of the inner 
spins are not fixed. The system is visualized in figure [|. 

The energy of the system as function of the two inner spins is H — —4 — 02 + 03 — <J\Oi. The 
energies for the four possible states are shown in the following table. 



02 


03 


H 


+1 


+1 


-5 


+1 


-1 


-5 


-1 


+1 


-1 


-1 


-1 


-5 



So the ground state of the system is threefold degenerated. It can be described by the rule: 
The inner spins are not fixed, but if 03 = +1 then 02 = +1 must hold too. So the ground 
state is formally given by 

V = {S, V X ,V 2 , T}, S = {1}, Vx = {2}, V 2 = {3}, T = {4} (8) 

and 

5 ={(2,1)} (9) 

The ground state can be displayed as a graph. The nodes represent the clusters and the 
edges represent the relation. Because for 'real' systems the clusters are too large, only the 
order parameter of each cluster is given in the graph, instead of enumerating all spins. For 
the inner clusters, where the order parameter can take two values, the value for the cluster 
orientation s(Vk) = —1 is displayed. So the example system is represented by the graph 
shown in figure ^. 

3 Calculation of ground states for all B, A 

Let's examine a Cluster C in a DAFF with magnetization mc = vn} c — ttiq (m^m^ > 
denote the number of spins in up and down direction). The cluster is supposed to exist as 
a coherent unit for an finite i3-interval [B\, B2\- Without loss of generality all other clusters 
can be supposed to keep their orientations inside this interval. 

The cluster is connected to its environment by Nq satisfied and Nq unsatisfied bonds. The 
part of the cluster energy describing its connection to its surface and to the external field is 

E C (B) = (N£ - Ng)J - m c B = (N£ - N^)J - (mj, - m l c )B (10) 

(NQ,NQ,m} c ,m\j depend itself on B.) For a ground state E C (B) — > min holds. Let's sup- 
pose that the cluster flips its orientation at B = B' J e [B±, B^, That implies that the 
magnetization of the cluster is nonzero, because the environment of the cluster is supposed 
to be unchanged, and otherwise there would not occur a jump. Since the number of satis- 
fied/unsatisfied bonds and the number of up/down spins exchange, when a cluster is flipped, 
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the energy of a nipped cluster changes its sign. That implies E C (B J ) = 0, because E C (B) < 
for all B 6 [B± , B2] ■ So the cluster can take both orientations at B = B^ . That leads to 

B>/J=^^f (11) 
m c — m c 

(NQ,NQ,m c ,mQ all from B < _B J or all from B > B 3 .) Because the cluster C is part of 
cluster S or T, these values are not directly accessible, we want to express this "jump-field" 
B 3 using characteristic values of the graphs R immediately "before" (B < B 3 ) and "after" 
(B > B 3 ) the jump. Because the cluster does not contribute to the degeneracy (which holds 
only for clusters with zero magnetization and Nq = N c ) it moves during the jump from the 
(super)cluster S to T or vice versa. So we get (ms,mT'. magnetizations of S,T) 



\<m c - m c \ = \m T {B < B 3 ) - m T {B > B 3 )\ (= \m s (B < B 3 ) - m s (B > B 3 )\) (12) 

Using Ng(B < B 3 )-N C (B < B') = N C (B < B^)-N^(B > B j ) = N U (B < B^)-N U (B > 
B^ ) (the rest of the system does not change, N u — total number of unsatisfied bonds in the 
system) and because B^ / J > we have as result: 

■ \N^B<Bi)-N"{B>Bi)\ 

1 \m T (B < Bi)-m T (B > &)\ { ' 

For a RFIM the formula describing the relevant terms of the energy of a cluster (similar to 
©)is 

E e ({B t }) = (N% - Nq)J - J2 (14) 

A given realization of a RFIM is determined by the local field values hi := -Bj(A = 1), so 
Bi = Abi. For the ± A- RFIM the bi is just the sign of the random field. Defining the total 
random field sign of the cluster by Ec := X)iec we § e t 

E C (A) = (JVg - Nq)J - E c s(C)A (15) 



Comparing with (10) we see that Ecs(C) plays the role of the magnetization and A is 
equivalent to the external field B of a DAFF. A RFIM cluster moves at a jump from the 
super-cluster S to T or vice versa, so we have |£c| = |£t(A < A J ) — St (A > A J )| and 
I s(C) I = 1. So the field A J where the cluster flips is calculated by: 

■ |JV"(A< AJ)-JV"(A> AQ1 

1 |Er(A < Ai) - S T (A > A^)| 1 1 

For illustration in Figures ^ to || the ground states of a sample DAFF system (L = 10, 
p = 0.5) are displayed for 9/7 < B < 1.6. In addition to the order parameter of each cluster 
(upper value) its magnetization (lower value) is given. 

A cluster with zero magnetization and equal number of satisfied and unsatisfied bonds can 
take both orientations over a range of external fields. For 9/7 < B < 1.5 these holds for 
clusters C, D. The degeneration of the system has degree 2x2 = 4, because each cluster C, D 
can take independently two orientations. 

At B^ = 1.5 a cluster J (an aggregation of the small clusters Ji, J2, J3) reverses its orientation. 
Here the number of total unsatisfied bonds is N u (9/7 < B < 1.5) = 39 and A^(1.5 < B < 
1.6) = 48 (not visible from the graphic representation, but easily calculated using the spins 
from clusters S, T, because all other clusters have — N%<). The magnetization of T 
changes from itit(B < B^) = 42 to itit(B > _B J ) = 48 (figure ||). Using formula ( [l3|) we get 
no lj — l 48-39 l =2 = 15 

D I J ~ |48-42| ~ 6 — L '° 

Each of the clusters G, Ji, J2 can take independently two orientations (—states). The clusters 
C, E alone can take three combinations of orientations because s(E) — +1 implies s(C) = +1. 
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The clusters Ji,D,F alone can take 4 states. Because s(E) = +1 implies s(Ji) = +1 we 
have for all 5 clusters together 4 + 4 + (4 — 1) = 11 states. So the system is 2 3 xll = 88 fold 
degenerated. 

For 1.5 < B < 1.6 the cluster J has been absorbed by the cluster S. Now the system has a 
degree of degeneracy of 2x3x3 = 18. Note that the environment of the cluster J has changed 
during the jump as well. But only clusters V with my = and N v = N v are affected. So 
the formula for B 3 /J is still correct. 

Using the formulas ( |l3| , |l6| ) the complete behavior of a DAFF/RFIM can be easily calculated. 
We show how the procedure works for a DAFF, the RFIM is similar. One starts with a 
certain value for B = Bq (There must no jump occur, otherwise the formulas are not valid) 
and calculates the graph R(Bq). Then B is increased in macroscopic steps AB up to the 
value B\ where the calculated graph R{B) differs from R(B ) the first time. Then the jump- 
field Bi can be calculated from the graphs at B , B\. One has to check if the jump is really 
the first jump occurring for B > Bq. This is done by calculating the graph R(B 3 _) for a field 
B 3 _ which is infinitesimal smaller than B 3 and comparing with R(Bq). If they are the same, 
the jump B 3 is really the first for B > B 3 . If not, a new jump-field B 3 is calculated using 
R{Bq) and R(B 3 _) and at B 3 may not be a jump at all. This procedure can be iterated until 
the jump next to Bq is found. 

It can be easily seen that this algorithm really produces the next jump starting from Bq. 
Suppose we have two jumps B\ < B 3 , between Bo and B\. Then B 3 is given by (for a DAFF) 

\m T (B ) - m T (Bi)| 

Because N u and mj are monotonic in B, the effect on that values of the both jumps sum 
up. So we have (B 3 + slightly larger than B 3 ) 

B'U = lNU{Bo) - jV "( g -)l + \ NU ( B +) ~ ggOl (lg) 
\mT{B )-m T (B 3 i)\ + \mT{B 3 + )-mT(B 1 )\ 

From B\ < B\ and because n\jd\ < n^ldi implies n\jd\ < (n\ + n%) /{di + efe) < "2/^2 we 
get 

\N U (B ) — N U (B 3 )| . \N u (Bi) - A^ n (Si)| 

B 3 J=± y —^- <B 3 /J < y -M y -Jdl = B 3 J 19 

\m T (B ) - m T (B 3 _)\ \m T (B 3 + ) - m^B^ 

By induction follows that the algorithm finds the next jump from Bq for all number of 
jumps between Bq and B\. The next jump from a given jump B 3 is found by restarting 
the algorithm at Bq := B 3 + . Beside this linear search a bisection method starting with 
Bq = B st art and B\ = B en d should be possible. But we implemented the linear search, 
because it is more straight forward. Taking into account, that the calculation of the graphs 
R requires B, J to be integers, the linear search of all states R in a range [B sta rt/ Jo, B enc i/ Jq] 
can be algorithmically displayed as follows: 
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algorithm all ground states {B start , B end , J ) 
begin 

B <— B s tart 

AB <- Jo 
stop <— false 
while ( not stop ) 
begin 

B ^B _ 
calculate R(B ) 
samc_R <— true 

while (same_R) /* Search for first state different from R(Bq) * / 
begin 

if(same_R AND B > B end ) then 

stop <— true; 
else 
begin 

B <- B + AB 
calculate i?(B) 
AB <- 2 * AB 
if( B(B ) ^R(B)) then 
same_R <— false 

end 
end 

if( stop) then 
continue 

found <— false 
Bi <- B 

while( not found) /* search for next jump */ 
begin 

& <_ ( int ) \ \ N y R a \- N " { %l * Jo + 0.5l - 1 

- v I y \m T (Bo)-mT(Bi)\ v J 

calculate R(B 3 _) 

if( R(BL) =R{B )) begin 

found <— true 
else 

B 1 ^Bi_ 

end 

There is a jump at B^ <- |iV u (B ) - iV u (Bi)|, J J <- |m T (B ) - m T (Bi)| 
B <- (int) [J * Bi/Jj + 0.5] + 1 
AB <— difference between last two jumps 
end 
end 



Technical remarks: 

• Jo should be chosen sufficiently large. For the DAFF and the ±A-RFIM J > N 
guarantees to find all jumps. Since the values for bi can take all real values for the 
Gaussian RFIM there is always a small probability of missing jumps. So we rounded 
the 6i-values to three significant digits and used Jo > lO 3 ^. 

• The field increment AB used for searching the next state which is different from R(Bo) 
is set to the difference between the last two jumps found (initially J ). During the 
search AB is doubled in each iteration, so jumps lying close to each other do not slow 
down the following calculations. 

• The states R(Bq),R(Bi),R(BL) must not be exactly at jumps in order to guarantees 
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that the algorithm works. This is confirmed by checking if there are inner clusters with 
nonzero magnetization. If this is true the field B is in-/decrease about 1 up to a non- 
jump field (not shown in the algorithmic representation). For the DAFF/ ±A-RFIM 
this is guaranteed to be done in one step. 

• For the DAFF /± A-RFIM the states R at the jumps (B j ,Jj) can be calculated. For 
the Gaussian-RFIM that is not always possible directly, because the values of Bi are 
rounded off. 

4 Results 

We performed calculations of ground states over large B. A-intervals (values are related to 
J = 1) for system sizes L = 8 to L = 48 for the DAFF (p = 0.5), ±A-RFIM and the 
Gaussian-RFIM. The number of samples per system type and size are shown the following 
table. 



L 


DAFF 


±A-RF1M 


Gaussian-RFIM 


B e [0, 61 


Ae [0,61 


A e [0,41 


8 


850 


450 


250 


12 


850 


450 


120 


16 


450 


250 


58 


24 


450 


177 


32 


32 


50 


43 


16 


48 


50 







For the L = 48 samples we reduced the i3-interval to [0.5,1.5]. For the L = 32 Gaussian- 
RFIM we reduced the A-intcrval to [1.0,2.8]. 
In figure || the antiferromagnetic order parameter per spin 

i 

of two DAFF systems with L = 8 and L = 16 over the whole range B E [0,6.5] is displayed. 
Because of the ground state degeneracy many states are possible. Here the two extreme 
values, the order parameter can take, are shown. For the small system for B < 1 two states 
are possible, because its total magnetization is zero. The larger system has this feature only 
for B < 0.5. For both systems the range of possible a- values is relatively small for 1 < B < 2 
and little bit larger for 2 < B < 4. The range is the largest for integer values / of B, because 
spins having neighbors with total magnetization I can take both orientations. The main 
differences between the two systems are: The larger system exhibits more jumps and the 
order parameter decreases more quicker than for L = 8. 

The discrete structure is clearly visible. In both systems the order parameter is high for 
small fields and goes stepwise to zero. 

Because of the step-structure it is difficult to obtain 'critical' magnetic fields by calculation of 
susceptibilities and then use finite size scaling to go to the L — > oo limit. Instead we proceed 
as follows, to get first informations about how the order parameter changes. We calculated 
for each sample the barrier field value B a<ao where the order parameter a falls the first time 
below ao, i.e. 

B a<aa = min{B\a(B) < a } (21) 

In figure [7] the average B a<ao is plotted against the inverse system size 1/L for ao = 
0.2, 0.5, 0.8. Just to get a first impression we extrapolated to L — > 00 by fitting straight 
lines B a<ao (l/L) = c/L + B a<ao {0)- The fit seems to be good, but may have no physical 
meaning. From the figure we can see that S a <o.5(0) = 0.90 ± 0.02, so the transition to 
the domain state occurs for B c < 1 ! This is in contradiction to previous results from MC 
calculations of systems of size 60x60x61 || where B c « 1.4 was found. Because we want 
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to rule out the influence of the fit function, we checked the result by calculating the ground 
states for ten L = 160 systems, at B = 0.9999 where we got an average order parameters of 
0.22 ± 0.11. So the infinite system has no long range order for B > 0.9 ! From figure we 
can't see, whether the [-B a <o.8j ^ a< o.2] interval remains finite or becomes infinitely small for 
the L — oo system. More calculations and larger systems are needed. 
In figure ^| diagrams are shown for the magnetization per spin 

i 

of two sample ±A-RFIM. 

The impression is similar to the DAFF, but the order parameter goes to zero at higher fields 
and there are more jumps in the system because there are more spins in a RFIM than in a 
DAFF of the same size L. Diagrams of Gaussian-RFIM look similar but there are even more 
jumps in a system of equal size and the decrease of the order parameter is less rapid. 
In analogy to the DAFF we calculated the barrier fields A m<mo for the random field systems. 
They are display in figure ^| 

We got the values A m<0 . 5 (0) = 2.35 ±0.01 (±A) and A m<0 . 5 (0) = 2.49 ±0.04 (Gaussian) for 
the infinite systems. The value for the ±A system is consistent with previous calculations 
§• 

For the Gaussian-RFIM the values are less reliable, because we could calculate only few 
systems. The Gaussian systems exhibit many jumps and the number of ground states which 
have to be calculated is about six times the number of jumps. So many ground states have 
to be calculate to analyze the system, resulting in more computer time needed. 
For all types of systems the number of jumps jfcj increases with system size. Figure ^ 
shows this dependence. The number of jumps increases roughly like — c * L b with 

b = 1.23 ± 0.04 (DAFF), b = 1.76 ± 0.01 (± A- RFIM), b = 2.91 ± 0.03 (Gaussian-RFIM, 
only jumps in A E [1.0, 2.8] are counted, the number of jumps in [0.0,6.0] is about four times 
higher !) This means that the behavior of the system gets smoother with increasing size, i.e. 
the infinite system looses the discrete structure on finite B, A-scales. 

We measured also the number of jumps in the B, A-interval where the order parameter falls 
from 0.8 to 0.2 for each system. Even this number increases although the interval is reduced 
with increasing system size. We performed a similar fit to the case above and obtained 
b = 0.90 ± 0.02 (DAFF), b = 1.40 ± 0.06 (±A-RFIM) and b = 2.18 ± 0.02 (Gaussian-RFIM). 
In the last part we focus on the degree of degeneration (dod) . In figure [n] the dod of the 
sample DAFF systems is displayed as function of the field B. The dod is not calculated 
exactly, but estimated very accurately^ by an algorithm found in Jl^| . 

The degree of degeneracy is there very large, where the external field takes integer values. 
For non integer values the dod takes the largest values for B £ (3,4). 

We found that there are only small connected subgraphs in R of small height (the maximum 
length path in the graph) at fields B where no jump occurs. So the spin-clusters contributing 
to the degeneracy are flipping almost independently, i.e. the ground state structure is quite 
simple. That means that the number of states are not reduced very much from 2#™, where 
#n is the number of inner clusters. 

We could find that the average dod can be approximated very good by a function of #n and 
of the number of edges #e in R over the full i?-range: 

dod = 2* n a* e (23) 

This means, that for #e = all inner clusters are independent, so the dod depends expo- 
nential on the number of nodes. The edges put constraints on the nodes, so the number of 
possible states is reduced. So it is possible to describe the ground state structure by a single 

2 For the graph types which occur here, the estimate seems to be exact: For small systems the value is the 
same. For large system the type of the graph is the same, as we will see. 
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value a. We measured a ~ 0.79 ± 0.03. This value is quite independent of the system size 
and the external field^. 

Figure [l^ shows the dod for the two sample ± A-RFIM systems. Here the dod takes its largest 
values at the integer-fields A = 3,4,6 and in (3,4). For the Gaussian-RFIM no degeneration 
occurs except the twofold degeneration at the jumps. This is to be expected, because it 
is very unlikely that the sum of the random fields Ec is zero for a cluster C . The dod of 
the ± A-RFIM can also be approximated very well by equation (p3|). We got an average 
a = 0.7582 ± 0.0002. 

Since the dod is mainly determined by the number #n of (inner) clusters and the number 
of edges, we calculated these values for the DAFF/± A-RFIM at B = A = 3.9999 for sizes 
up to L = 1 60 (from 2560 realizations for L = 12 to 10 realizations for the largest systems). 
The result is shown in figure 

For both system types all numbers have a L 300 dependence on L. The number of edges is 
about 3 to four times smaller than the number of clusters. For the DAFF about 4.6 percent 
of the spins are contained in the inner clusters, for the ±A-RFIM this value is only 1.8 
percent. These values are almost independent of the system size. Because the number of 
clusters/state depend linearly on the system size, this shows again that the ground state 
of the DAFF/RFIM consist of many independent small clusters. So the number of clusters 
increases, but not the size of the clusters and not the degree of connection between the 
clusters. This means that the degree of degeneracy depends exponential on the size of the 
system. 



5 Conclusions 

In this paper we have presented an algorithm for calculating exact ground states of random 
field systems and diluted antiferromagnets in polynomial time for all fields B, A. The com- 
plete degeneracy of all states is representable by a relation of (anti-)ferromagnetic clusters. 
The use of new fast algorithm allowed us to calculate systems of sizes up to 160 3 . 
We found that the transition to the domain states occurs for the DAFF at B c < 1 in con- 
tradiction to previous results. The result for the ±A-RFIM A c = 2.35 ± 0.01 was consistent 
with previous calculations. For the Gaussian RFIM we found A c = 2.49 ± 0.04. 
The order parameters of the systems are stepwise constant functions of the field. The number 
of jumps inside a given B, A interval increases faster than L, indicating that the infinite 
system has a rather smooth behavior. 

The ground state degeneracy of the DAFF and the ±A-RFIM is very high, especially when 
the field is integer valued. Apart from these values the degree of degeneracy is at highest in 
the domain state. The degree of degeneracy increases exponentially with system size. The 
degeneration is rather simple, because the ground states consist of many small, only loosely 
coupled clusters. Not more than 5 resp. 2 percent of all spins contribute to the degeneracy. 
The Gaussian-RFIM is not degenerated at all, except at the jump fields, where one cluster 
can take two orientations. 

Using the methods presented and referenced in this paper many more information about 
the ground state structure of DAFFs and RFIMs will be obtained in the future. A detailed 
analysis of the numbers, sizes and forms of the clusters forming the ground states and type of 
connections between them is possible. Because the systems get smoother with increasing size, 
finite size scaling techniques should be applicable, to gain more insight into the transition 
from the long range order to the domain state. 
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Captions 

1. Example RFIM system with four spins. 

2. Ground state of example RFIM given by the relation R. 

3. Ground state R of sample DAFF (L = 10,p = 0.5) for 9/7 < B < 1.5. 

4. Ground state R of sample DAFF (L = 10, p = 0.5) for B = 1.5. 

5. Ground state R of sample DAFF (L = 10, p = 0.5) for 1.5 < B < 1.6. 

6. Antifcrromagnetic order parameter of sample DAFFs (L = 8, p = 0.5) and (L = 16, 
p = 0.5.) for < B < 6.5. 
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7. Barrier field B a<ao (l/L) for L = 8, . . . , 48 (a = 0.2,0.5,0.8) and fit functions B a<ao (l/L) 
c/L + B a<ao {0). 

8. Antifcrromagnetic order parameter of sample RFIMs (L = 8) and (L = 16) for < 
A < 6.5. 

9. Barrier field A m<mo (l/L) for L = 8, . . . , 32 (m = 0.2,0.5,0.8) and fit functions 

(1/L) = c/L + (0). 

10. Average number of jumps in B e [1,2.8] (Gaussian-RFIM) resp. [0,6] (±A-RFIM, 
DAFF) as function of system size in a double logarithmic scale and fit functions 
#j(L) = cL b . 

11. Degree of degeneracy of sample DAFFs (L = 8, p = 0.5) and (L = 16, p = 0.5) for 
< B < 6.5. 

12. Degree of degeneracy of sample RFIMs (L = 8) and (L = 16) for < A < 6.5. 

13. Number of clusters #n, edges #e as function of system size L for the DAFF and the 
±A-RFIM together with fits of form f(L) = cL b . 
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